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Fluxes of Atmospheric Neutrinos and Related Cosmic Rays 

T.K. Gaisser a * 

a Bartol Research Institute, University of Delaware, 
Newark, DE 19716, USA 

The atmospheric neutrino beam simultaneously spans a range of pathlengths from ten to ten thousand kilome- 
ters, which correspond respectively to downward- and upward-going neutrinos. As with any neutrino oscillation 
experiment, also in this case the interpretation of the data depends on a detailed knowledge of the neutrino beam. 
The ingredients are the primary spectrum of cosmic-ray nucleons, the geomagnetic fields in which the charged 
particles propagate and the properties of interactions of hadrons in the atmosphere. In this talk I review the 
status of calculations in light of the recent evidence for neutrino oscillations from Super-Kamiokande jl[ . 
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1. Introduction 

When cosmic ray protons and nuclei interact in 
the atmosphere, the secondary cascades include 
neutrinos from decay of pions, muons and kaons. 
Production of these neutrinos depends on the lo- 
cal zenith angle because of the competition be- 
tween decay and interaction of the parent mesons 
in the tenuous atmosphere. A simple geometric 
construction |2) shows that a trajectory from be- 
low with nadir angle 9 has the same zenith angle 

9 on the other side of the earth. Therefore, since 
neutrinos with E <C 10 5 GeV are virtually unat- 
tenuated by the earth, the flux of atmospheric 
neutrinos would be up-down symmetric in the ab- 
sence of neutrino oscillations except to the extent 
that the isotropy of the primary cosmic rays is 
distorted by the geomagnetic field. Variation of 
the neutrino flux with azimuth is a consequence 
only of the geomagnetic field (the "East- West" ef- 
fect) even in the presence of oscillations (because 
within a given band of zenith angle the distri- 
butions of neutrino pathlengths and energies are 
independent of azimuth) . This fact allows || an 
important check of the systematics of the Super- 
K analysis which I discuss in §2. 

Typical altitudes of production of the neutri- 
nos are between 10 and 20 kilometers, so the 
distribution of neutrino pathlengths ranges from 

10 km for vertically downward neutrinos (neutri- 



nos that originate directly overhead) to ~ 10 4 km 
for upward neutrinos from below. The range of 
neutrino energies for contained or partially con- 
tained events is from sub-GeV to multi-Gev. For 
neutrino-induced upward, throughgoing muons it 
extends to ~ 1000 GeV. Thus the atmospheric 
neutrinos have a range of pathlength over neu- 
trino energy 1 < L km / E GeV < 10 5 . 

The 7r — > /i — > e decay chain is the predominant 
mode of production of atmospheric neutrinos in 
The sub-GeV to multi-GeV range. This leads to 
the basic prediction of 



v u + v u 



(1) 



for E, y < 1 GeV. The ratio decreases as energy 
increases because muons are increasingly likely 
to reach the surface before decaying. Compari- 
son of decay length to energy-loss length in the 
earth leads to the conclusion that virtually all 
muons that reach the ground stop before decay 
(or capture) occurs. Therefore muons that reach 
the ground do not contribute the neutrinos with 
energy high enough to contribute even to the sub- 
GeV sample (p e > 100 MeV/c). 

In contrast with the expectation of Eq. |l|, sev- 
eral experiments find 



R = 



(p - like/e- Uke) data 
(jj, — like/e — like)MC 



0.65, 



(2) 
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which is equivalent to (y e + v^jiv^ + i>^) w 0.77. 
Ingredients that enter the denominator of the ra- 
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tio of ratios in Eq. || include the calculated neu- 
trino flux, the cross sections for neutrinos to in- 
teract in the quasi-elastic and various multi-prong 
channels and the detection and reconstruction ef- 
ficiencies in the detector. The subject of this talk 
is the neutrino fluxes. 

Three independent calculations [p|-pr| have 
been compared and analyzed Q in order to iden- 
tify and evaluate the sources of uncertainty in 
our knowledge of the atmospheric neutrino flux. 
Here I concentrate on comparison of the calcu- 
lations of Honda et al. 10 13] with the "Bar- 



tol fluxes" 0,|T^1 , because these two have been 
used by the experimental groups P,p|Jl^] for anal- 
ysis of their data. The calculations of Refs. |l4|] 
are extensions of the calculations of 
respectively to high (> 3 GeV) and to 
low (< 200 MeV) energy. In addition, a more 
realistic treatment of the geomagnetic cutoffs is 




used 1 17 1, which reduces the calculated neutrino 
fluxes in the sub-GeV range by about 10 per cent. 
Honda et al. jl3| also extended their calculation 
to high energy. At present the Bartol neutrino 
fluxes and the Honda et al. fluxes in the GeV 
range agree within 5% in magnitude as well as 
ratio. This level of agreement in magnitude is, 
however, smaller than the systematic uncertain- 
ties, as I will discuss in Sections 3 and 4. 

In sections 3 and 4 I discuss the primary fluxes 
and the treatment of hadronic interactions, both 
of which influence the spectrum and shape of the 
neutrino spectra. In the Super-K analysis, the 
overall normalization is treated as a free parame- 
ter because of the large uncertainty in the normal- 
ization of the primary spectrum. There are re- 
cent measurements of the primary spectrum that 
should in principle allow one to reduce this source 
of uncertainty. 

As emphasized by Perkins Jli| , muon fluxes 
high in the atmosphere arc directly related to the 
neutrino fluxes, being produced by the same pri- 
mary spectra and from the same interaction pro- 
cesses. In §5 I discuss how measurements of the 
flux of muons high in the atmosphere are being 
used to check the overall normalization and shape 
of the closely related neutrino flux. In the con- 
clusion I list the various approximations common 
to the present calculations and how they might 



be expected to affect the results. 

2. Geomagnetic effects 

Propagation of a cosmic-ray nucleus through 
the geomagnetic field depends only on its gyrora- 
dius and hence on the magnetic rigidity, 



R 



A x pc/(Ze) 



(3) 



Here A and Z are the mass and charge of a 
nucleus of momentum-per-nucleon p. Low en- 
ergy particles at low geomagnetic latitudes can- 
not reach the atmosphere to produce secondaries. 
Since energy per nucleon is the important quan- 
tity for production of secondaries, nuclei become 
relatively more important compared to protons 
at low geomagnetic latitudes because of the fac- 
tor A/Z w 2 in Eq. |. 

Both neutrino flux calculations 13 



1, 



use 



geomagnetic cutoffs obtained by the standard 
method of backtracking antiprotons through the 
geomagnetic field to determine the cutoffs for a 
particular location. For example, in the calcu- 
lation of Ref. [0, which is used in Ref. |T^| , 
antiparticles are injected at 20 km altitude on 
an outward trajectory. If the trajectory reaches 
30i?e before it travels 500i? ffi and without in- 
tersecting the surface of the earth, then it is as- 
sumed that positive particles of the same rigidity 
can reach the atmosphere from that direction. 

For the location of Super-K we have compared 
the cutoffs used in Ref. |l3| with those of Ref. ]17y 
used for the calculation of Refs. |l4|,[l5|]. The cut- 
off maps are very similar, but with some notice- 
able differences toward the east, where the cutoffs 
are slightly higher in Ref. |l7j . 

At low geomagnetic latitudes such as Kamioka, 
average cutoffs are higher locally (i.e. for cos- 
mic rays entering the atmosphere above the de- 
tector) than for the opposite hemisphere (i.e. for 
the cosmic rays entering the atmosphere on the 
other side of the earth, which give rise to upward- 
going events). The opposite is the case for a de- 
tector at a high geomagnetic latitude, such as 
Soudan. There the local cutoffs are negligible 
in the sense that essentially all cosmic-rays from 
above with sufficient energy to produce pions and 
contribute to the flux of neutrinos can reach the 
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atmosphere to interact. Upward events originate 
from the atmosphere over the entire hemisphere 
below each detector. Since the average over a full 
hemisphere is similar from any viewpoint, the up- 
ward/downward ratio should be greater than one 
at Kamioka but less than one at Soudan. 

Fig. H illustrates the situation. The pair of 
curves labelled (A) shows the distribution of pri- 
mary cosmic-ray energies that would contribute 
to the sub-GeV signal in Super-Kamiokande if 
there were no geomagnetic cutoff at all. The solid 
curve is for solar minimum and the dotted one for 
solar maximum. The middle pair (B) is the cor- 
responding response from below, which would be 
similar if Super-K were moved to Soudan. The 
rightmost pair of curves (C) is the response for 
downward sub-GeV events at Super-Kamiokande. 
What is plotted is proportional to the event rate 
per logarithmic interval of primary energy, so in 
each case the area is proportional to the signal. 
Thus the upward/downward ratio at Super-K is 
B/C > 1. If Super-K were located at Soudan, the 
ratio would instead be B/A < 1. (The method 
used to simulate "sub-GeV" events is described 
in Ref. §.) 
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Figure 1. Response for sub-GeV muon-like events 
in Super-K. 



As neutrino energy increases the up-down 
asymmetry from the geomagnetic effect dimin- 
ishes. For this reason, the Super-Kamiokande 
group have emphasized the multi-GeV event sam- 
ple in their search for neutrino oscillations Jl|,|l9| . 
On the other hand, the full data set contributes 
to the evidence for oscillations. Therefore it is im- 
portant to note H that the geomagnetic effects 
themselves provide a way of testing the integrity 
of the entire analysis chain that is independent of 
whether or not there are oscillations. 

At the low geomagnetic latitude of Kamioka 
there is a pronounced east-west effect on the cos- 
mic radiation. Cutoffs are significantly lower for 
positive particles from the west than from the 
east. For example, the trajectory of a 20 GeV 
antiproton injected toward the east from above 
Super-K at 70° from the zenith would be bent 
down by the geomagnetic field and intersect the 
surface of the Earth, while the same antiproton 
injected toward the west would escape from the 
geomagnetic field. In other words, the cutoff for 
protons with zenith angle 70° from the east at 
Super-K is > 20 GeV. For directions closer to 
the horizon the cutoff from the east approaches 
50 GV. In contrast, for directions above the hori- 
zon from the west the cutoff is 5 to 10 GV at 
Kamioka. 

The excess of primary cosmic rays from the 
west at Kamioka produces a corresponding east- 
west asymmetry of the low-energy neutrino flux 
and hence of the sub-GeV event rate. There is 
a much smaller, but still non-negligible asymme- 
try for the multi-GeV event sample || . Since the 
east-west effect is an azimuthal asymmetry, it is 
independent of oscillations; oscillation effects de- 
pend on neutrino pathlengths, which vary with 
zenith angle but are independent of azimuth. 

Figure |2Q| compares the azimuthal depen- 
dence of the Super-K data (0.4 < pi ep ton < 
2.0 GeV/c, single ring events in 22.5 kton fidu- 
cial volume) with expectation. The solid line uses 
the neutrino fluxes of Ref. fill a nd the dashed 
line the calculation of Ref. [|l4| , |l5| . Although the 
fits are equally good (y^/d.o.f. ~ 1 for all four 
comparisons) , the geomagnetic effect is somewhat 
more pronounced with the Bartol neutrino flux 
than with the flux of Honda et al. |||. 
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We have made some diagnostic tests to investi- 
gate the source of this difference and a similar dif- 
ference between the two calculations that shows 
up in the zenith angle dependence of sub-GeV 
events. The difference arises in part from the dif- 
ference in cutoffs mentioned above, but also from 
the difference in primary spectrum, as discussed 
below. 
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Figure 3. Summary of primary spectra. 



Figure 2. Azimuthal dependence: calculated and 
observed by Super-K. 



3. Primary spectrum 

Both the normalization and the shape of the as- 
sumed primary spectrum have important conse- 
quences for the calculation of the neutrino fluxes. 
The normalization propagates directly through to 
the event rate. The assumed spectral index af- 
fects the shape of the neutrino energy spectrum 
in an obvious way, but it also affects the angular 
dependence through its interaction with the ge- 
omagnetic effects. Thus, a softer spectrum will 
lead to more pronounced geomagnetic effects be- 
cause a larger fraction of the event rate comes 
from lower energy primaries, which are most af- 
fected by the geomagnetic field. 



Fig. |^ is a summary of measurements of spectra 
of protons, helium and the CNO group of nuclei, 
compared with the primary spectra of Honda et 
al. Q (solid lines) and the spectra used in Ref. 
(dashed lines). From Fig. |l|, it is apparent 
that 5 < E < 50 GeV/nucleon is the most im- 
portant region of the primary spectrum for sub- 
GeV events. The harder spectrum of Ref. |L3| in 
this energy region, coupled with the geomagnetic 
field, contributes significantly to the fact that the 
geomagnetic effects are somewhat smaller in the 
neutrino spectra of Ref. |R| than in Ref. |Q , as 
mentioned in the previous section. 

A marked feature of the plot for hydrogen is 
the fact that the data of Webber [glj (shown by 
the open circles in Fig. |^) are significantly higher 
than those from the LEAP experiment |^2j and 
other more recent experiments in the same energy 
region. The difference is outside the error bars, 
indicating a systematic effect. The most recent 
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result is from the BESS detector |23|, and other 
recent experiments (Refs. |24|,|2f|) are included in 
the BESS compilation. Generally (with the pos- 
sible exception of the measurement of Ref. Pq| ), 
the interpretation of which is complicated by an 
unusually strong level of solar modulation) all the 
recent experiments are consistent with the LEAP 
results. In fits to their data the Super-K group 
have treated the overall normalization of their 
rates as a free parameter. The primary spectra 
and their potential consequences for interpreta- 
tion of the data on neutrino interactions are dis- 
cussed more fully in talks given at the Satellite 
Symposium 

4. Yields 

It is important to note that the primary spec- 
trum is not the only source of uncertainty in 
the normalization and shape of the energy spec- 
trum of atmospheric neutrinos. Uncertainties in 
the yields of pions and kaons in interactions of 
hadrons with nuclei of the atmosphere are also 
important. Not all of phase space is covered in ac- 
celerator measurements with nuclear targets. For 
sub-GeV events the important range of beam en- 
ergies is from a few GeV to several tens of GeV 
(see Fig. [l]). In this energy range the atmo- 
spheric cascades are dominated by interactions of 
nucleons, and nearly all neutrinos are from the 
7r — > /i — > e decay chain. 

Existing measurements with beam energies 
around 20 GeV and light nuclear targets measure 
pions only above 3 or 4 GeV ^9j,|o) , and there are 
significant differences in how the lower energy pi- 
ons are represented in the different neutrino flux 
calculations, as discussed in Ref. p2| . The pion 
multiplicities, and the momentum distributions 
as reflected by the spectrum-weighted moments 
for pion production, are highest in the calculation 
of ||,[l||,[l5| ■ This compensates to some extent for 
the higher assumed proton spectrum of Ref. [fl3| 
with the result that the calculated neutrino fluxes 
(comparing Refs. |Q and Jll|) differ by less than 
either the primary spectrum or the yields. 

Yields in a new calculation of Battistoni et 
al. a re intermediate between those of Refs. [[)) 
and . 



5. Muons 

The same primary spectra and the same 
hadronic interactions determine both muon and 
neutrino fluxes. Therefore, comparison with mea- 
surements of muons high in the atmosphere offers 
a way to check directly the neutrino fluxes. The 
most important range of altitudes for pion decay 
is 10 to 25 kilometers, which corresponds to at- 
mospheric depths of ~ 20 to ~ 200 g/cm 2 . 

Many of the same detectors referred to above in 
connection with recent measurements of the pri- 
mary spectrum have also been used to measure 
the muon spectrum during ascent through the at- 
mosphere and on the ground. The calculations of 
Refs. JlJ] (and jl3|]) compare reasonably well with 
the measurements of the MASS experiment [^2| , 
although there is a relative excess of muons be- 
low 1 GeV in the calculation. On the other hand, 



a recent comparison between |14 and the HEAT 
measurements of muons J33| showed better agree- 
ment in the shape of the spectrum but with an 
overall excess of the calculation relative to the 
data of as much as 50% in some bins. 

Measurements on the ground and at float alti- 
tude necessarily have better statistics than data 
obtained during ascent. It is possible that some 
of the discrepancies referred to above could be 
a consequence of the short exposures during as- 
cent. Both for MASS §§ and HEAT @ there is 
a tendency for better agreement between calcula- 
tion and measurement at float and at the ground 
than during ascent. This is an active area with 
further potential for reducing uncertainties in the 
flux of atmospheric neutrinos. 

There are interesting possibilities with the 
muon measurements for probing details of the cal- 
culations. For one thing, muon fluxes at float 
altitude reflect directly the primary spectrum 
and the properties of pion production in single 
nuclcon-nitrogen interactions with no intervening 
cascading. 

A more interesting possibility arises from the 
fact that in some cases the same detector has 
been exposed at different locations with differ- 
ent geomagnetic cutoffs. The MASS experiment 
has been flown both in Northern Canada (essen- 
tially no cutoff) and from Ft. Sumner, NM where 
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the vertical cutoff is « 5 GV. The BESS detec- 
tor has measured the muon charge ratio on the 
ground in northern Canada and in Japan. The 
low-energy behavior of the ratio is quite different 
in the two locations. Below ~ 1 GeV the 
ratio decreases toward 1 in Japan, which can be 
understood as a consequence of the high local ge- 
omagnetic cutoff. There are two effects: 
First, with a high cutoff, heavy primaries are rela- 
tively more important because they have a higher 
rigidity for a given energy per nucleon. Pro- 
tons produce more positive than negative pions 
(and hence more and vice versa for neutrons. 
Thus, enhancing the contribution from nuclei, 
which carry the neutrons, suppresses the muon 
charge ratio slightly. Secondly, vertical /i + at the 
ground have followed trajectories from slightly 
east of vertical (where the cutoffs are higher), 
whereas vertical fi~ will have come slightly from 
the west where more of the primary spectrum 
reaches the atmosphere to produce secondaries. 
This enhances the negative relative to the posi- 
tive muons and hence reduces the fi + j ' pT ratio 
preferentially at low energy where the bending is 
more significant. 

6. Conclusion 

Present calculations |ll|,[l4| include several ap- 
proximations: 

• They are one-dimensional; i.e. all neutrinos 
are assumed to follow the direction of the 
primary nucleon that produced them. This 
approximation has two effects: 

f. There should be some loss of parti- 
cles that are produced at large an- 
gle. Given the momentum involved, 
as compared with the typical trans- 
verse momentum of produced pions, it 
is straightforward to check that this ef- 
fect should be small for neutrino events 
in Super-K. 

2. Bending of charged particles in the 
atmosphere is not followed. This is 
perhaps the most important effect to 
check Q because it is systematic. 
As explained above, the vertical muon 



charge ratio is reduced when the cut- 
offs are high. There is a correspond- 
ing decrease in the v e jv e ratio (and 
an increase in that part of the v^jv^ 
ratio that comes from muon decay). 
Because ov > cr&, the calculated ra- 
tio of electron-like to muon-like events 
will decrease with respect to the one- 
dimensional calculation. This correc- 
tion will therefore make the anomaly 
of Eq. (||) somewhat more pro- 
nounced. 

• The superposition model has been used in 
Rcf. (TJ] for interactions of nuclei. Within 
the framework of a standard multiple scat- 
tering picture, this approximation can be 
shown to give a good account of the distri- 
bution of first interactions of each nucleon. 
ft will, however, lead to some overestimate 
of the multiplicity of pions in the target 
fragmentation region. 

• The cascades are propagated to sea-level all 
over the globe. In particular, the exact ter- 
rain over the detector (i.e. the mountain 
in the case of Super-K) has been neglected 

This is negligible for muon neutri- 
nos from pions, which decay high in the 
atmosphere. From the pathlength distribu- 
tions of Ref. [El it is possible to estimate 
the size of the effect of this approximation. 
To take an extreme case of a 4 km overbur- 
den (^NUSEX), the neutrinos from muon 
decay overhead are overestimated by about 
10%, leading to a ~ 5% overestimate of the 
calculated (y e + t? e )/(v^ + v^) ratio. 

• The calculations are based on parametriza- 
tions of data in limited regions of phase 
space. Interpolations and extrapolations 
introduce some level of uncertainty. The 
yields of Ref. | {l4| a re at the high end of a 
spectrum, with Jl l| the lowest and |l3| and 
PH in between. 

At least two groups are embarking on three- 
dimensional calculations. The Italian group ]3l[ | 
has published a short account of their plan with a 
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comparison of their one-dimensional results with 
those shown in Ref. |§. They use FLUKA j|§ 
with various hadronic interaction models at dif- 
ferent energies. The authors of Refs. || and 
Jli), together with Coutu, are also pursuing this 
goal. Although effects are generally expected to 
be small, in view of the importance of the ex- 
perimental results, a greater level of detail in the 
calculations is warranted. 

Acknowledgements. I am grateful to Ed 
Kearns for providing me with Fig. 2 [ ^o| and 
to Todor Stanev for reading the manuscript and 
for collaboration on this work. I thank M. Gold- 
haber and M. Spiro for useful conversations, and 
M. Honda for very helpful exchanges of informa- 
tion about the calculations of Refs. |^0| , |l3| . 

REFERENCES 

1. Y. Fukuda et al, Phys. Rev. Letters 81 (1998) 
1562. 

2. D. Ayres et al, Phys. Rev. D29 (1984) 902. 

3. Paolo Lipari, Todor Stanev & T.K. Gaisser, 
Phys. Rev. D (to be published) 



ph/9803093 



astro- 



Y. Fukuda et a l Phys. Lett. B 433 (1998) 9 
(sub-GeV) and |hep-ex/980500l (to be pub- 
lished) (multi-GeV). 

5. T.K. Gaisser & Todor Stanev, Phys. Rev. D57 
(1998) 1977. 

6. R. Becker-Szendy et al. (1MB Collaboration), 
Phys. Rev. D46 (1992) 3720 and references 
therein. 

7. Y. Fukuda al. (Kamiokande Collaboration) 
Phys. Lett. B335 (1994) 237 and references 
therein. 

8. W.W.M. Allison et al. (Soudan Collabora- 
tion), Phys. Lett. B391 (1997) 491. 

9. Giles Barr, T.K. Gaisser & Todor Stanev, 
Phys. Rev. D39 (1989) 3532. 

10. M. Honda, K. Kasahara, K. Hidaka & S. Mi- 
dorikawa, Phys. Lett. B248 (1990) 193. 

11. E.V. Bugaev & V.A. Naumov, Phys. Lett. 
B232 (1989) 391. 

12. T.K. Gaisser, M. Honda, K. Kasahara, H. 
Lee, S. Midorikawa, V. Naumov & Todor 
Stanev, Phys. Rev. D54 (1996) 5578. 



14 



15 



16 



17 

18 

19 
20 



21 



22 
23 



25 
26 
27 

28 

29 
30 

31 

32 



13. M. Honda, T. Kajita, K. Kasahara & S. Mi- 
dorikawa, Phys. Rev. D52 (1995) 4985. 
Vivek Agrawal, T.K. Gaisser, Paolo Lipari & 
Todor Stanev, Phys. Rev. D53 (1996) 1314. 
T.K. Gaisser & Todor Stanev, Proc. 24th Int. 
Cosmic Ray Conf. (Rome) vol. 1 (1995) 694. 
M. A mbrosio et al (MACRO Collabora- 
tion), hcp-cx/9807005] (throughgoi ng muons) 
and M. Spurio, hep-ex/9808001 (stopping 
muons). 

Paolo Lipari & Todor Stanev, Proc. 24th Int. 
Cosmic Ray Conf. (Rome) vol. 1 (1995) 516. 

D. H. Perkins, Astroparticle Physics 2 (1994) 
249. 

T. Kajita, to appear in Proc. of Neutrino98. 
C. McGrew (Super- Kamiokande Collabora- 
tion), to appear in Proc. Int. Conf. on High 
Energy Physics (Vancouver, 1998) and Super- 
K paper on azimuthal dependence (forthcom- 
ing). 

W.R. Webber, R.L. Golden & S.A. Stephens, 
Proc. 20th Int. Cosmic Ray Conf. (Moscow) 
vol. 1 (1987) 325. 

E. S. Seo et al. (LEAP) Ap.J. 378 (1991) 763. 
S. Orito, T. Sanuki et al. (BESS), to appear 
in proceedings of the conference " New Era in 
Neutrino Physics" , Tokyo, June 11/12 (1998). 

24. W.R. Menn et al. (IMAX) Proc. 25th Int. 
Cosmic Ray Conf. (Durban) vol. 3 (1997) 409. 
M. Boezio et al. (CAPRICE), Ap. J. (to be 
published) . 

W.R. Webber et al. (MASS), Ap.J. 380 (1991) 
230. 

M. Honda, to appear in proceedings of the 
conference "New Era in Neutrino Physics", 
Tokyo, June 11/12 (1998). 
T.K. Gaisser, to appear in proceedings of the 
conference "New Era in Neutrino Physics", 
Tokyo, June 11/12 (1998). 
T. Eichtcn et al, Nucl. Phys. B44 (1972) 333. 
J.V. Allaby et al, CERN Yellow Report No. 
70-12 (unpublished). 

G. Battistoni et al, Proc. 5th TAUP Conf., 
Gran Sasso (1997). 

M. Circella, C.N. De Marzo, T.K. Gaisser & 
Todor Stanev, Proc. 25th Int. Cosmic Ray 
Conf. (Durban), vol. 7 (1997) 117. 
33. Stephane Coutu (HEAT Collaboration) to 



8 



appear in Proc. Int. Conf. on High Energy 
Physics (Vancouver, 1998). 

34. Maurice Goldhaber has emphasized this point 
(private communication) . 

35. I thank M. Spiro for raising this question. 

36. A. Fasso, A. Ferrari, J. Ranft and P.R. Sala, 
Proc. of the 3rd Workshop on Simulating Ac- 
celerator Radiation Environment, SARE-3, 
KEK-Tsukuba, May 7-9, 1997, (H. Hirayama, 
ed.) KEK Report 97-5, p. 32 (1997). 



